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Abstract: 

Within the QCD-improved parton model and assuming the factorization theorem 
to hold in the helicity basis and for higher twist contributions, we show how non 
zero single spin asymmetries in hadron-hadron high energy and moderately large px 
inclusive processes can be obtained, even in massless perturbative QCD, provided the 
quark intrinsic motion is taken into account. A simple model is constructed which 
reproduces the main features of the data on the single spin asymmetry observed in 
inclusive pion production in pp collisions. 



1 Introduction 



Single spin asymmetries in high energy and moderately large pp inclusive hadronic 
processes have recently received much attention, both experimentally [fj]-|4j and the- 
oretically P|-||13||. Whereas they are expected to vanish at leading twist in massless 
perturbative QCD, higher twist effects might still be important in the kinematical 
region of the available data and may give origin to non zero values. Among the 
attempted explanations, quark-gluon correlations || j|, [L3| and transverse k± effects 
in the quark distribution || |7J or fragmentation functions JTTJ have been considered. 

We analyse here the large single spin asymmetries observed in the collision of 
transversely polarized protons off unpolarized protons, with the production of pions 
with pt values up to 2 GeV/c, 

p^p -> nX . (1) 

Several experimental results are available for such processes and show clear 

patterns in the dependence of the spin asymmetry 

dcr^ — da^ 

An = da^ + 6*i (2) 

on xp = 2pL/yfs, where pp is the pion longitudinal momentum in the pp cm. 
system, and yfs is the total cm. energy. The proton spin is perpendicular to 
the scattering plane and there exist data for ir^ and 7r°; proton spin orientations 
parallel to the scattering plane would give, by parity invariance, zero single spin 
asymmetries. A striking dependence of A N on p T which was reported at xp ~ [|3| 
is now believed, after a more careful analysis of the data |14[|, not to be observed 



anymore; A^ixp ~ 0) shows no sign of dependence on pp and remains consistent 
with zero in the whole range of pp explored so far, (0 <pp ^4) GeV/c |T|| . 

Some previous data fli~5]-||17|| seem to be in disagreement with the findings of Refs. 



I]- 1|, [|Tj]; in particular, they do not show sizeable values of A N for n |T6|, [L7j], 
show strong pp dependences for 7r° |TJj and 7r + [17 and have contradictory results 



on xp dependences [|16| |17j. However, all these data come from experiments at much 
smaller energy (proton beams of 13.3, 18.5 or 24 GeV/c) than those of Refs. -[[§], 
|L4| (200 GeV/c) and we do not consider them here because they might be outside 
the range of applicability of perturbative QCD. 

The usual QCD description of large pp inclusive production is based on the fac- 
torization theorem, according to which the cross-section for the production, say, of a 
large pp pion in the collision of two protons is given by the convolution of an elemen- 
tary cross-section with the number densities of quarks and gluons inside the protons 
(distribution function) and the number density of pions inside a quark or gluon (frag- 
mentation function). The elementary cross-section contains all the dynamical and 
quantum-mechanical information on the constituent interaction, whereas the dis- 
tribution and fragmentation functions are phenomenological ways of modeling the 
non-perturbative long distance physics: they can be measured in other processes 
and their large Q 2 evolution is given by perturbative QCD. 
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We know that the above hard scattering scheme works well for unpolarized pro- 
cesses and indeed it has been tested in many experiments. It has also been general- 
ized to the polarized case jl8|l , so that it may be applied to the description of several 



processes involving polarized hadrons fI5 |; however, the existing spin data do not 
allow yet a definite test of its validity. 



In Section 2 we adapt the formalism of Refs. || and [IB to the case of the single 
spin asymmetry (0). In order not to obtain a zero result higher twist effects have 
to be introduced; this can be done at different stages and several suggestions or 
attempts have been proposed in the literature HH HI- Single spin effects 

in the elementary reactions alone are bound to be proportional to a a m q j^fs WA - 
where m q is the quark mass, and are then expected to be negligible at high energies, 
even taking constituent quark masses into account |j, [2(| . Spin effects might 



then be present in the distribution or fragmentation functions: the former has been 
suggested by Sivers || [7| and the latter by Collins |J, whose idea has been further 
developed and applied in Ref. [[□]]. Qiu and Sterman [§] have used both higher order 
elementary interactions and higher twist distribution functions to predict a sizeable 
single spin asymmetry in large pt direct photon inclusive production, p^p — * 7X. 

The approach we describe here is equivalent, although derived in a different way, 
to the suggestion of Refs. f| [7|] and supports it; we then discuss (Section 3) a simple 
model which implements the idea and gives very good agreement with the data. The 
approach of Refs. [|, [7j has been criticized in Ref. || on the ground of violating 
the time reversal invariance of QCD. This is true only at leading twist order, if soft 
initial state interactions between the colliding protons are neglected, which need not 
be the case in Sivers or our model; we will further comment on this in Section 2. A 
short conclusion is given in Section 4. 

2 The single spin asymmetry in the hard scattering scheme 

Let us then consider the process (1), supposing the initial protons moving along 
the z-axis and choosing xz as the production plane; the incoming proton is polarized 
parallel (f ) or opposite (|) the y-direction so that, in the helicity basis, 

lt> = ^(|+>+<|-» (3) 

li> = ^(l+>-<|-»- (4) 

According to the QCD factorization theorem the differential cross-section for 
the hard scattering of a polarized proton with spin j (and similarly for spin \) on 
an unpolarized target proton, resulting in the inclusive production of a pion with 
energy E n and three-momentum p n , p^p — > nX, can be written as |9|, [18|, |i~9 



E n da pl ! p^ wX 1 r 1 

dx a dx b - (5) 



d3 P7t 2 a,b,c,d \ a A' a ;X b ;\ c ,\' c ;\ d J Z 
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where f a / p t (x a ) is the number density of partons a with momentum fraction x a inside 

/ T 

the polarized proton [similarly for fb/ P {xb)] and p a x x , (x a ) is the helicity density 
matrix of parton a inside the polarized proton p\ The Ma c ,a ;A a ,A 's are the helicity 
amplitudes for the elementary process ab — > cd; if one wishes to consider higher 
order (in a s ) contributions also elementary processes involving more partons should 
be included. D^jf c (z) is the product of fragmentation amplitudes 

where the X stays for a spin sum and phase space integration of the undetected 
Tx,x x 

particles, considered as a system X. The usual unpolarized fragmentation function 
D n / C (z), i.e the density number of pions resulting from the fragmentation of an 
unpolarized parton c and carrying a fraction z of its momentum is given by 

D„c{z) = l -Y. D l7>^)- (7) 

K 

For simplicity of notations we have not indicated in Eq. (|5|) the Q 2 scale dependences 
in / and D; the variable z is related to x a and Xf, by the usual imposition of energy 
momentum conservation in the elementary 2 — > 2 process |23|; we have skipped, for 
the moment, some spin independent kinematical factors^, but we explicitely kept 
the factor 1/2 to remind that an average has been taken over the helicities of the 
unpolarized parton b (quark or gluon). 

Eq. (|5p holds at leading twist and large pt values of the produced pion; the 
intrinsic k± of the partons have been integrated over and collinear configurations 
dominate both the distribution functions and the fragmentation processes; one can 
then see that, in this case, there cannot be any single spin asymmetry. In fact, 
total angular momentum conservation in the (forward) fragmentation process [see 
Eq. ([])] implies A c = A^,; this, in turns, together with helicity conservation in the 
elementary processes, implies A a = A^. If we further notice that, by parity invariance, 
D w J c c does not depend on A c and that J2x b ,x c ,x d |^A c ,A d ;A a ,Aj 2 is independent of A a we 

remain with J2x P\ P \ = 1- Moreover, in the absence of intrinsic k± and initial state 
interactions, the parton density numbers j a j v \ (x a ) cannot depend on the proton spin 
and any spin dependence disappears from Eq. (|), so that 

da prp ~™ x - do- plp ^ x = . (8) 

Eq. (|5|) can be generalized with the inclusion of intrinsic k± || and this can 
avoid the above conclusion; for example 0, the observation of a non zero k± of a 

1 Equation (Q) holds as an equality if the elementary amplitudes are normalized so that 
(1/4) Ea .x b .x .x d \Mx a ,x b ,x a ,xf = (lM(da/di) 

a ■ b c ' d ° a 
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final particle C with respect to the axis of the jet generated by parton c does not 
imply any more A c = X' c and allows a non zero value of the asymmetry 



(9) 



The above asymmetry (^|) is related to the so called Collins [11], |24|, [25| or sheared 
jet [26] effect; it requires the measurement of the azimuthal angle of the outgoing 
hadron around the jet axis, but, apart from a small sin0 dependence, it is a leading 
twist effect and it depends on some non perturbative quark fragmentation analysing 
power. When integrating over the azimuthal angle the effect might not entirely 
disappear because of some dependence in the elementary parton interaction. This 
idea was exploited in Ref. |TT) where, essentially, the parton c is produced in the 
forward direction and the final hadron px is due to its transverse k± inside the jet. 
One cannot expect such a model to work at large px- 

Another possible k± effect, suggested by Sivers || |7j, may originate in the dis- 
tribution functions. To see how this comes out from the general scheme we rewrite 
Eq. (^) taking into account the parton intrinsic momentum in the number density 

fa/p- 



1 

2 



E 



E / d 2 kj_ a dx a dx b - 



(10) 



P\ a P \> a fa/pi(x a , k ±a ) f b/p (x b ) Mx c ,x d -,x a ,x b M*,^.^ ZT^ C (z) , 

where / denotes the k± dependent number density. 

We can now argue, as in the previous case when no k± was taken into account, 
that angular momentum, helicity and parity conservation eliminate all dependences 
on the parton helicities in Eq. (0); however, a dependence on the hadron spin 
may remain in f a / p i(x a , k± a ), analogously to the Collins effect in the fragmentation 
process. Then one has |6|, [/fl: 

E n da p]p ^ nX E n da pLp ^ nX 



A„,A' 



1 

4 

x 



d 3 p n 

E E 



d 3 p n 
d k i f, doCn doch 



(ii) 



fa/p^{ x a, k±a) fa/pi-\ x aik-La) fb/p{%b) 



M; 



A„,A d ;A £t ,A b 



D 



tt/c( 



Several comments are now in order. 

There is a new quantity which appears in Eq. (|TT|): 



/pt 



k± a ) 



ai k±a) 



fa,X a /pt {x a , k± a/ 



E [fa,K/P^ X 

,X a /p^{Xaik_\_a) fa,X a /pl{Xa) 



(12) 
(13) 
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where f aX / p f(i) (x a , k± a ) is the number density of partons a with helicity A a , momen- 
tum fraction x a and intrinsic transverse momentum k± a in a transversely polarized 
proton [with spin f or J, according to Eq. (§) or (§)]. Eq. ([13]) follows from Eq. 
(|12|) by rotational invariance and explicitely shows that A f a / p r(x,k±) = when 

fcj_ = 0. 

This new quantity can be regarded as a single spin asymmetry or analysing power 
for the p 1 ' — > a + X process; if we define the polarized number densities in terms of 
distribution amplitudes as 

fa,xj P r(x a ,k ±a ) =T \Q a ^ XaA {x a ,k ±a )\ 2 (14) 

J Xp,\ Xp 

then we have, in the helicity basis, 

A N f a/p T(x a , k± a ) = T J2 2lm K P ,A ■+(Xa,k ±a ) G%?* X ._(x a ,k ±a ) 

J x p ,x Xp Xa p p 

= 2I a + /p (x a ,k ±a ). (15) 



Eq. ( Jl5| ) simply follows from Eqs. ([12]) and (|T4| ) via Eqs. (|3]) and (0J) and shows 

a? &_La)- 



the non diagonal nature, in the helicity indices, of I+^(x a , k 



Collins H has argued that a non zero value of I^/l(x, k±) is forbidden by the time 
reversal invariance of QCD; his argument is based on the analysis of the non diagonal 
matrix elements of the leading twist quark operator ipj+ip, whose diagonal matrix 
elements are related to the distribution functions f q / P (x, k±). For such operator his 
argument is correct and indeed time reversal invariance forces the matrix elements 
between proton states with different helicities to be zero. In a more physical language 
this amounts to say that single spin asymmetries for the process p 1 ' — > qX are 
forbidden by parity and time reversal invariance, which is true. However, as we said, 
initial state interactions (like soft gluon exchanges) between the incoming protons 
must certainly occur, and, at least in cases in which their neglect gives a zero result, 
one should consider them and relate the distribution functions to the inclusive cross- 
section for the process P1P2 — > qX; single (transverse) spin asymmetries are then 
certainly allowed (as the problem we are studying here confirms) via time reversal 
invariant scalar quantities like ^^vpaVxV^ s \ ■ 

These soft gluon and initial state interactions which correlate partons from dif- 
ferent hadrons and allow a non zero value of I+£ can only survive at higher twist; it 
has been shown [[Il| that at leading twist-2 the proof of the factorization theorem for 



the unpolarized case - with the cancellation of all soft gluon contributions - holds 
for the polarized well. Eq. (|TT]) , as it will be shown in the next Section, only 

gives higher twist contributions; assuming a non zero value of 1+^ in the factorized 
structure of Eq. ( |TT]) amounts to assume the validity of the factorization theorem 
beyond leading twist. 

In the operator language this approach has been advocated by Qiu and Sterman 
I! who use generalized factorization theorems valid at higher twist and relate non 
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zero single spin asymmetries in pp collisions to the expectation value of a higher 
twist operator, a twist-3 parton distribution, which explicitly involves correlations 
between the two protons and combines quark fields with a gluonic field strength. 
However, they still consider only collinear partonic configurations so that, in order 
to obtain non zero results, they have to take into account the contributions of higher 
order elementary interactions. 

Our function I+l(x a , k± a ) introduced in Eqs. ( |15D and (fT2|) can be considered 
as a new phenomenological quantity which takes into account the non perturbative 
long distance physics, including initial state interactions, and plays, for single spin 
asymmetries in pp collisions, the same role plaid by the distributions functions 
fa/pi x a) in unpolarized processes. This function is zero in the absence of parton 
intrinsic motion, but, for k±_ ^ 0, allows non zero single spin asymmetries even 
taking into account only lowest order perturbative QCD interactions among the 
constituents. Measurements of the asymmetries supply information on 1+1, like 
measurements of unpolarized cross-sections supply information on f a / p . 

In the next Section we introduce a simple parametrization of 2+^, based on our 
knowledge of the distribution functions, and show that it can reproduce with good 
accuracy the data on the single spin asymmetry 



da p]p ^ nX - da plp ~* nX da p]p ^ nX - da plp ^ nX 

Afj = = 

da p]p ^ x + da pLp -* wX " 2 da unp 



(16) 



3 A simple phenomenological model 

Inserting Eqs. (|TT| ) and (]T5f) , with the proper kinematical factors [23|, into Eq. 
dT6|) yields 

An _ T,a,b,c,dJ ' dx a dx b d 2 k ±a I+-{x a , k±g) f b / P {x b ) [da/di(k ±a )] D n/c (z)/z 
Ea,b,c,dl dx a dx b f a /p(x a ) f b/p (x b ) (da / di) D w/c (z)/z 

where da/dt is the unpolarized cross-section for the elementary constituent process 
ah — > cd; notice that the dependence on k± a has to be kept into account in such a 
quantity, otherwise the numerator of Eq. (|17D would vanish upon integration over 
k± a due to the fact that I+I(x a , k± a ) is an odd function of k± a , see Eq. (fTBD . The 
constituent momentum fraction z carried by the pion can be expressed in terms of 
x a and x b by requiring energy-momentum conservation in the elementary scattering. 

All quantities appearing in Eq. ( |T?D are either theoretically or experimentally 
known, with the exception of the new quantity I+£(x a ,k± a ), which, in principle, 
is measurable via the single spin asymmetry A^. However, in order to give an 
estimate of and see if we can obtain reasonable values within our approach, we 
parametrize here 1+^ in a most simple way. 

We assume that the dependence of 1+1 on k± a is sharply peaked around an 
average value k^_ a = (fc^) 1 / 2 , value which may depend on x a ; the x a dependence of 
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which does not originate from the k± a dependence is taken to be of the simple 

form 

N a x a a a {l-x a f\ 

so that we approximately have 



Jd 2 k ±a I^(x a ,k ±a )^(k ±a ) 



(k ±a ) x >o 
k° 



M 



a N a x^(l-x a 





La) 




da 


(+k± a ) 


da 


di 


di 


da 




da . 
~ ~dl^ 



C 1_J 



(18) 



(19) 



where M is a hadronic mass scale, M~l GeV/c. A numerical estimate of k\jM 



can be found in Ref. [^7] and can be accurately reproduced by the expression 



n '±a 
M 



0.47o£ 68 (l 



X ( 



\0.48 



(20) 



which we adopt in our calculations. Such value of k\ a enters in the computation 
of [da /dt (k± a ) — da/dt (—k± a )]; notice that such a quantity is 0{k^/px) @, hence 
the contribution of Eq. |T3 is a higher twist one. 

In order to give numerical estimates of the asymmetry dl7f) we still need explicit 
expressions of the unpolarized distribution functions, f a ,b/ P , and the fragmentation 
functions, D^/ c . At this stage we have only considered contributions from u and d 
quarks inside the polarized proton, which certainly dominate at large xf values, that 
is a = u, d in the numerator of Eq. (|TTD. Instead, we have considered all possible 
constituents in the unpolarized protons, with k± = 0, and all possible constituent 
fragmentation functions. We have taken f q ,q >g / p from Ref. p^| , D^/g^ from Ref. 



29] and D n / g from Ref. |Kj • Given the very limited pt range of the data we have 
neglected the QCD Q 2 dependence of the distribution and fragmentation functions. 
By using Eqs. ( |19"D and ( |20"D into Eq. §L7\), together with the unpolarized f a ,b/ P 



and D^/c functions, we remain with an expression of A n still dependent on a set of 



6 free parameters, namely N a ,a a and f3 a (a = u,d), defined in Eq. (|18|j . We have 
obtained a best fit to the data 0, shown in Fig. 1, with the following values of the 
parameters: 



u 5.19 2.79 4.15 



d -2.29 2.11 4.70 



(21) 
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As the experimental data [TJ] cover a px range between 0.7 and 2.0 GeV/c we have 
computed An at a fixed value px = 1.5 GeV/c. Our asymmetry decreases with 
increasing pt and increases at smaller pp; however, we do not expect our approach 
to be valid at Pt smaller than, say, 1 GeV/c. We will further comment on the Pt 
range of our computation in the conclusions. 

Notice that the above values ( PH ) are very reasonable indeed; actually, apart 
from an overall normalization constant, they might even have been approximately 
guessed. The exponents a u ^ and (3 u> d are not far from the very naive values one can 
obtain by assuming, as somehow suggested by Eqs. (p~5|) and that I^(x a ) 



-/p,+( x a)i where f a ^-y Pt+ (x a ) denotes, as usual, the number density 
of quarks with the same (opposite) helicity as the parent proton. Also the relative 
sign and strength of the normalization constants N u and turn out not to be 
surprising if one assumes that there might be a correlation between the number of 
quarks at a fixed value of k± a and their polarization: remember that, according to 
SU(6), inside a proton polarized along the y direction, P y = 1, one has for valence 
quarks P£ = 2/3 and P* = -1/3. 

It might appear surprising to have approximately opposite values for the 7r + and 
7r~ asymmetries, as the data indicate, and a large positive value for the 7r°; one 
might rather expect ~ for a 7r°. However, this can easily be understood from 
Eq. ( |l7|) which we simply rewrite, for a pion 7r*, as A l N = Af t /T>' t , if one remembers 
that from isospin symmetry one has: 

AHVc=^(ArVc + Ar-/ C ) • (22) 

Eq. (flTD show that the relation ( p2]) also holds for J\f° and X>°, so that 



1 + 



A N = JV +JV = A N • (23) 

i + v + 

It is then clear that A N ~ — A N implies A N ~ only if T>~ ~ V + , i.e. if the 
unpolarized cross- sect ions for the production of a 7r~ and a tt + are approximately 
equal. This is true only at xp ~ 0. At large xp the minimum value of x a kine- 
matically allowed increases and the dominant contribution to the production of 7r + 
and 7r~ comes respectively from f u / p (x a ) and fd/ P (% a ) [ see the denominator of Eq. 
dl7|)1. It is known that fd/ P (x a )/ f u / p (x a ) — > when x a — > 1; this implies that 
T>~ /T> + decreases with increasing xp, so that at large xp we have T>~ /T> + 1 and 
— — 2T>~ /T> + ) ~ A^f. Such a trend emerges both from the experimental 

data and our computations. 

4 Conclusions 

We have applied the QCD hard scattering approach, based on the factorization 
theorem in the helicity basis, to the description of single spin asymmetries in inclu- 
sive production of pions, with p T values in the 1 to 2 GeV/c range, in the scattering 
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of polarized protons off unpolarized ones. Such region of px is a delicate one; even if 
we assume, somewhat optimistically, that the hard scattering scheme is suitable for 
the description of these processes, we expect that, due to the moderate pr values, 
higher twist contributions may still be sizeable and important. Indeed, the leading 
twist contribution to the single spin asymmetries is zero and the large experimental 
data should be explained via non leading terms. 

In evaluating them one has to assume that the factorization theorem still holds 
at higher twist level; moreover, one can only take into account a few out of the 
many higher twist contributions and corrections to the simple hard scattering for- 
mulae. Here, we have considered the role of the parton intrinsic k± in the initial 
polarized proton 0, |7j, while neglecting other possible effects due, e.g., to transverse 
k± in the fragmentation process. In this our model can only be regarded as a phe- 
nomenological approach to the description of the otherwise mysterious large single 
spin asymmetries. Further application of the same model should test its validity. 

We have shown how single spin asymmetries can be different from zero and 
originate from lowest order perturbative QCD interactions, provided some non per- 
turbative and intrinsic k± effects are properly taken into account. In the helicity 
basis, suitable for the use of the factorization theorem, this non perturbative long 
distance information is contained in the function I+I(x,k±), which has a simple 
parton model interpretation in the transverse spin basis, Eqs. fl!2| ) -(|T5| ) . 

Detailed data on single spin asymmetries in pp inclusive interactions would yield 
information on I+£(x, k±), similar to the information gathered on the unpolarized 
structure functions; a knowledge of I+H from some process could then be used to 
predict spin effects in other processes. Here we have shown how a simple and realistic 
parametrization of 1+1 can easily explain the data on single spin asymmetries for 
the process p*p — > ttX at large energy and px — 2 GeV/c. Our model clearly 
exhibits the observed increase with xf, at fixed px values, of the magnitude of the 
asymmetries. A more detailed application of our approach to other processes, like 
p^p — > 7rX, p^p — > 7X and Trp"< — > ttX is in progress this should help in 

assessing the relevance and importance of our estimates. 
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Figure caption 

Fig. 1 Fit of the data on An p|, with the parameters given in Eq. (|2l|); the upper, 
middle and lower sets of data and curves refer respectively to 7r + ,7r° and ir~ . 
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